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Abstract 

In this work we consider the dynamical Casimir effect for a massless scalar field - under Dirichlet 
boundary conditions - between two concentric spherical shells. We obtain a general expression for 
the average number of particle creation, for an arbitrary law of radial motion of the spherical shells, 
using two distinct methods: by computing the density operator of the system and by calculating 
the Bogoliubov coefficients. We apply our general expression to breathing modes: when only one 
of the shells oscillates and when both shells oscillate in or out of phase. Since our results were 
obtained in the framework of the perturbation theory, under resonant breathing modes they are 
restricted to a short-time approximation. We also analyze the number of particle production and 
compare it with the results for the case of plane geometry. 
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I. INTRODUCTION 

The static Casimir effect, theoretically predicted in 1948 [1], consists in the attraction 
of two perfectly conducting, parallel plates, due to the distortion of the electromagnetic 
vacuum state. As a consequence of quantum fluctuations under the very presence of external 
constraints, such a distortion of the vacuum state emerges from the fundamental concepts of 
quantum field theory |2j . Casimir himself was the first to discuss the importance of spherical 
geometry in the distortion of the vacuum state [3|, proposing in 1953 a semiclassical model 
for the stability of the electron. In this model, the electron was assumed to be a perfectly 
conducting spherical shell carrying a total charge e, with the Coulomb repulsion balanced by 
an attractive Casimir force. However, as elucidated by Boyer [4| in 1968, the Casimir force 
in spherical configuration is repulsive, invalidating the attempt to explain the stability of the 
electron through the Casimir force. Nonetheless, the electron model proposed by Casimir 
provides evidence that one cannot predict whether the Casimir force will be attractive or 
repulsive before the whole calculation is carried out. Moreover, the development of the 
bag model of the hadrons in the early 1970s — describing hadrons as quarks and antiquarks 
confined inside of a spherical cavity — also stimulated the investigation of the Casimir effect 
with spherical geometry {5]. 

Beyond static geometries, the dynamical Casimir effect (DCE), arising from movable 
external constraints, probes even more deeply into the complexity of the vacuum structure. 
The dynamics of the geometry gives rise to a time-dependent Casimir force along with a 
dissipative component 61]. The mechanical energy dissipated by this 'vacuum viscosity' 
induces the most striking manifestation of the DCE, i.e., the mechanism of particle creation 
and annihilation. G. T. Moore was the first to proceed to the quantization of the radiation 



field in a cavity with movable perfectly reflecting boundaries [Tjl and, still in the 1970s, the 

HQ 

creation of photons from the nonuniform motion of the boundaries was predicted [8J, |9J. 
liater, it was also noticed that a sudden change in the refractive index of the medium 
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12j could also extract photons from the vacuum radiation field. 



Regarding the spherical geometry — the subject of the present work — it was proposed 
by Schwinger 13j that the DCE could provide a possible explanation for the sonolumines- 
cence phenomenon, discovered in the 1930s 14J. However, in spite of interesting results 
following Schwinger's hypothesis 15], the theoretical description of sonoluminescence still 
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remains controversial . We finally mention that the DCE with spherical geometry bears some 
similarity to the problem of particle creation in the expanding Universe, where the spatial 
portion of the metric is a hyperspherical surface with time-dependent radius [l6, 17, fla. 

In the present study we deal with the DCE for a massless scalar field confined between 
two concentric spherical moving shells, and present a general expression — for any law of 
radial motion of the shells — to compute the average number of particle creation. We note 
that the particular case of a unique spherical shell moving with a specific law of motion was 
studied in Ref. j^ol ]. After deducing an effective Hamiltonian for the problem, we compute 
the average number of particle creation by two distinct methods: by considering the time 
evolution of the density operator of the cavity field and also by computing the well-known 
Bogoliubov coefficients. Assuming, then, an oscillatory radial motion for the spherical shells, 
our results are applied to four different cases: when only (a) the inner or (b) the outer shell 
oscillates, apart from those where both shells oscillate (c) in phase or (d) out of phase. 

The present paper is organized as follows: in section II we present the field quantization 
and derive an effective Hamiltonian; in section III we calculate the average number of particle 
production, for an arbitrary law of motion; in section IV we specialize for four breathing 
modes and analyze our results; finally, in section V the present our concluding remarks. 



II. FIELD QUANTIZATION BETWEEN TWO MOVING CONCENTRIC 
SPHERICAL SHELLS 

To quantize a massless scalar field between the two spherical shells, we start from the 
action 



S 



dtd 3 x £(x) 



X - j dtd 3 x V0-^> 2 



where the Lagrangian density L enables us to evaluate the canonical momentum 

dC 1 , 



7r(r,t) 



d 



(v,t). 



(1) 



(2) 



By minimizing the action ([I]), we obtain the Klein-Gordon field equation 



(r,t)=0, 



(3) 
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where the cavity field 0(r, t) is subject to the Dirichlet boundary conditions 



(f)(r=n, 9, <p, t) = <f>(r = r Q , 9, if, t) = 0, (4) 

on the inner and outer spherical shells, with radii Ti and r D , respectively. The spherical 
geometry of the cavity leads us, naturally, to seek solutions for the cavity field and its 
canonical momentum in the form of spherical harmonic expansions 

oo I oo / — — 

<KM) = E E ^2cJ^F ls (r)[c lms (t)Y lm (e,ifi) + H.c], (5a) 

1=0 m =-l s=l ' ls 



oo / oo 



C V 2 

1=0 m=—l s=l 



7T(r, t) = -i E E - V ^T^M [Cb«(*)nn(e, ¥>) - ^-C] • (5b) 



By substituting the above expansions into Eqs. (j3J) and (Tj0), we obtain the differential 
equation 



r 2 dr \ dr 
under the boundary conditions 



I d^dF llW \ + ^_/0_M) |Fi<(rHOi (6) 



r 



= ri) = F is (r = r ) = 0. 

Moreover, as the radial functions are solutions to a boundary value problem, they automat- 
ically satisfy the orthonormality relations 

F ls (r)F ls ,(r)r 2 dr = 5 SjS ,. (7) 

As the solution of Eq. ([6]) is given by a linear combination of spherical Bessel functions 
of the first (ji) and second kind (nj), the boundary condition on the inner shell leads to the 
relation 

F ls {r) = Ms [ji {—) th {—) ~ Ji (—J th (—J , 
whereas that on the outer shell results in the transcendental equation 

31 (—) »i (—J - 31 (—J "1 (—J = °- ( 8 ) 

The index s in — which assumes discrete values and are not necessarily equally spaced 
- indicates the sth root of Eq. (151) . We also note that the derivation of the solution of 
the problem of the moving shells with dynamical boundary conditions follows directly from 



the replacement of the static r i(o ) by the dynamical radii r^ ){t), since all time dependence 
in the system arises from them. In Fig. 1 we have constructed a map of the solutions of 
the Eq. (jSJ) for some values of the numbers I and s. As we can see, for the case / = 0, 
the frequencies are equidistant, which does not occur for the case I 7^ 0. However, if we 
have r Q (t) » r Q (t) — r* (t) (i. e., when the radii of the shells are much larger than the 
separation between them), the solutions for all values of / approach to the solution for 1 = 0, 
i.e., ujis — > uj 0s - 



A. Canonical field quantization 



The canonical quantization of the scalar field <fi in Eq. (I5a| is performed — by promoting 
the coefficients c\ ms and c* lms to operators ai ms and a\, m , sl — through the construction of a 
field operator satisfying Eqs. ©-(jl]) and the equal-time commutation relation 



4>(r,t),n(r',t) 



(r,i),0(r',t)J = [7r(r,t),7r(r',t)] =0, 

where fr is the momentum operator associated with 7T. The above relations between the field 
operators automatically imply the bosonic commutation relations for the standard creation 
and annihilation operators 



a>ims(t), <4 mV (t) 



[aims{t),a Vm , s ,{t)] = aj ms (t),al, m , s ,(t) 



0. 



Through the time derivative of the quantum version of Eqs. ([5]), together with the equations 
for the cavity field ([3]) and its canonical momentum (T5]), we obtain the Heisenberg equation 
for the annihilation operators 



0>lr, 



5 (t) = -m s (t)ai ms (t) + } j Hi[s S '](t)aim S '{t) + Hi(ss')(t)<4(_ m)f> ,(t), 



(10) 



where /i z(ss /) (t) = [^w(t) + fii s > s (t)] /2 and fi l[ss/] (t) = [^ ss '(t) - Hi s > s (t)] /2 are the symmet- 
ric and antisymmetric parts, respectively, of the coefficient 



H>lss'{t) = lS ^\ $ss> + (1 - 8 SS >) 
ZUl s [t) 



r 2 F ls ,(r;t)F ls (r;t)dr. 



n{t) 



From Eq. fflOl) we directly obtain a 



■J 



Ims' 



B. An effective Hamiltonian 

Following the reasoning exposed in Ref. [21[, we next derive an effective Hamiltonian 
governing the evolution of the creation and annihilation operators, as given by Eq. ffTUj) . To 
this end, we consider the most general form of a quadratic Hamiltonian 

Heff = \fll^mm'ss'(t) a \ ms a \'m's' + fll'mm' ss'it) a \ms a Vm' V 

1,1' m,m' s,s' 

+fllWnm'ss'(f) a lmsa\i m i s i + fu' mm 'ss' (t) a lmsO-l'm's' , 

which governs the evolution di ms = (i/h) [H e ff, aims]- Comparing the evolution equations for 
aimsit) and a\ ms (t) obtained through the Heisenberg equation of motion with those following 
from Eq. (fTOl) . we obtain the effective Hamiltonian H e ff(t) = Ho(t) +V(t), where 

H (t) = H <*>,,(*) (a\ ms a lms + ~ J , (11a) 

l,m,s ' 

V(t) = MsAt) [ a lms' + <4(~m)s') a Ls ~ a lms (a«-m)a' + 4 m s') ■ ( llb ) 

l,m s,s' 

In what follows we compute the average number of particles created in a selected mode 
via two distinct methods: the density operator and the Bogoliubov coefficients. 



III. AVERAGE NUMBER OF PARTICLE CREATION 



A. The density operator 

From our effective Hamiltonian H e ff(t) we obtain, in the interaction picture, the density 
operator of the cavity modes 

Pit) = P(0) + {-Ji) J Q dt i I dt 2 • ■ • J Q dt n [Vjih), [Vj(t 2 ), ■ ■ ■ [Vj(t n ),p(0)]}} 
where 

Vi(t)=i^Y,Y,Piss>(t) 

l,m s,s' 

with hmsif) = a/ ms exp (-ifta(t)), a] ms (t) = a] ms exp (iQi s (t)), and Qi s (t) = /*dti ui 8 (ti). 

To compute the average number of particles created in a particular mode labeled by the 
quantum numbers (I, m, s), given by Ni ms (t) = Tr p{t)a\ ms ai ms , we go up to the second-order 
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aims' (t) + 4_ m)s ,(t)J aimsit - aimsit) (ai^s'it) + a\ ms ,{t) 



approximation in the velocity of the cavity walls, fj, ro -C c. Starting from the vacuum state 
p(0) = |{0}) ({0}|, we thus obtain 



N lms (t) 



1 /"* Z"* 1 

- / dii / dt 2 V VV/ii w '(ii)/i W (t 2 

<rr) ' <n-)" n n' n n' 



m' ,m" q,q' p,p' 



x a im"qai(-m")q>a\ ms a lms a\^_ m , )p ,a\ m , p |{0}) 

x exp [-i (fij p (ti) + Qip'ih) - ni q (t 2 ) - n v (t 2 ))] 

+ ({0}| ^m' P a/(-m')p'«L S a ^ a f(-m'')g' a L" g K }) 

x exp [z (fifp(ti) + ftip'ih) - VL lq {t 2 ) - VL w {t 2 ))}} . 



(12) 



Using the result 



<{0}| a i(-m')p' a lm'p IW) = ^ 

m" ,m^m' ,m^q' ,p'^ s,p^ si,qi ^m" ,m& —m! ,vnflq' ,pO s,p'^s,q 
3—m",m3m',m3q,p'3s,p3s,q' —m" ,m& —m! ,m^q,p^ s,p'Os,q' > 



the expression (}I2"]) reduces to the compact form 



(13) 



Ni ms (t) = ^ / dti^( s / s )(ti) exp {z [Oj s /(tl) + 
which, like the energy of a given mode Wj s , is the same for any value of m. 

B. The Bogoliubov coefficients 

In this section we c omp ute the average number of particle creation by means of the 
Bogoliubov coefficients [22j, defined as 



air, 



q=l 



q=l 



4m S (*) = ^ a *« g (*) a Lff(°) + As 9 (*) a i(-»™)</(0), 



(14a) 
(14b) 
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with the initial conditions ai sq (0) = 5 S)q and (3i sq (0) = 0. In what follows, we obtain the time 
derivatives aj ms (t) and a/ ms (t) directly from Eq. (fllj) : 

= ate g (*)aimfl(0) + ^ A s9 (04(-m) 9 ( )' ( 15a ) 

*L(*) = X) d ^(*)°U(°) +E^( t ) a '(-^(°)> ( 15b ) 

and also by substituting ai ms (t) and a| ms (t) from Eq. (|T4|) into Eq. (J9]), we get 

ai ms (t) = [-^ s (5 SjS 'Q;; S(? (t) + Hl[ ss ']Oils'q(t) + fl>l(sa') Pis' q{t)] a lmq(0) 

s',q 

+ [-iuJi s 5 s ^f3i sq {t) + fii [ss '](3w q {t) + ii Kssl) al, q (t)] af ( _ m)9 (0), (16a) 

s',q 
s',q 

+ [iui s 8 StS ,ffi sq (t) + LH\ss'}Pis> q {t) + M(ss>)Oci s >q(t)] a/ ( _ m)? (0). (16b) 

By equating terms in the expressions for cti ms (t) and aj ms (t) in Eqs. (fT5|) and ffTBT) we obtain, 
after some algebraic manipulation, two coupled differential equations for the Bogoliubov 
coefficients 

ais S '{t) = -iuji s (t)ai ss '(t) + ^2 [iii[ sq ](t)ai qs *(t) + Hi{sq){t)$ qs >{t)\ , (17a) 

q 

$iss'(t) = -iUl s (t)f3l ss >(t) Y [P>l[sq]{t)(3lqs>{t) + M(sq) (*)<**„,/ (*)] • (17b) 

Now, we expand these coefficients in powers of the coupling strength /^ ss /, such that 

oo 

«U0=e^ (t) j>W(f) (18a) 

A=0 

oo 

/M*)=e^ (t) (18b) 



A=0 
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where the factor e is introduced for convenience. By substituting Eq. (|18|) into Eq. 

(117p we derive the recurrence relations 



(A) 

a 



\(t) = [ d tl Yl \*M<h) e ~ tQlq(tl) a t> l) ^) + 

(19a) 

= T dti e^^i £ L N] (*i) e"^ < s 7 1} (ti) + ftw^) e*W*> ag^'ft)' 

(19b) 

When the initial conditions, given in the zeroth-order terms a[%{t) = S s>s > and /3^j/(t) = 0, 
are substituted back into Eq. (119bl) . we finally obtain the first-order solution 



/©(*)= / dhe'to-toM'Winwih). (20) 



We next compute the average number of particle creation from the expression Ni ms (t) = 
({0}| o\ ms (t)ai ms {t) |{0}) = ^2 S , \f3i tS ,s'{t)\ 2 , where |{0}) indicates the initial vacuum state of 
the cavity and t is the time interval during which the shells have been in motion. Up to 
second order in the coupling coefficients fii ss >, we do obtain the result 

N lms (t)^J2\risl(t)\ 2 ■ (21) 

s' 

The substitution of Eq. ( 1201) into Eq. ( 1211) finally gives exactly the same Eq. ( fl3i) . 



IV. PARTICLE CREATION BETWEEN HARMONICALLY OSCILLATING 
SHELLS 

In this section, we assume that the shells perform small harmonic oscillations described 

by 

r a {t) — f a [1 + e a sin (tut)] , a = i,o. (22) 

where e a < 1 and w stands for the frequency associated with the oscillating shells. By 
substituting Eq. (1221) into Eq. ( fl3l) we obtain, up to second order in e a , the result 
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E 



exp [i (ui ss/ + w)t]-l exp [i (ui ss > - w) t] - 1 x 



(u>l s3 > + w) (uJlss' — ZO) 

2 

X 



Yl C ?(ss>f a taZU , (23) 



where we have defined 



c iss> = 7; — ttt — ^ o Sfl / + (1 - <W) 4/ 7— / rF Zs /(r;0) — dr 



2a; Zs (0) dr a \ ^(0) J r <9r a 



+ (1 _ Wl/ ^W\ W )f^d r , (24) 

and ti^ss/ = cji s (0) + u>i s /(0). From Eq. (|23|) . we observe the occurrence of resonances when 
w = trigs', and for a given mode I, s, the average particle creation in the s'th resonance is 
given by 

lim N lms (t) ~ cf (ssl) r a e a tut (25) 

exhibiting a quadratic increase with time. 

For I = the transcendental equation (jHJ), has the analytical solution 

S7TC 

^(*)=^oi(*) = ro(t) _ ri(t) , 

relating the radii to the instantaneous field frequencies for mode s, implying the resonance 
condition w = (s + s') ^oi(O). The coefficients thus reduces to 



r * — ( 1 ~' 

c 0(ss') - I X J c 0(ss') - r _ r . ' 

and the average number of particle creation is given by 

ss' \ / e r - (-1) e^r 



Vss* 1 



lim iV Qms (i) ~ ( - - _ ) I v ^ ^ I (tut) 2 . (26) 



[s + s') 2 J \ r -ri 



Notice that since the lower bound of r Q — is \r e \ + |rj£j|, the maximum value of the second 
factor on the RHS of Eq. f[2"6"j) is 1. This is in agreement with the fact that the Casimir effect 
is more pronounced at small distances between the shells. Moreover, the effective velocities 



zue r — (— l) s+s we^i = v a — (— l) s+s Vi play an important ro 
process, which is compatible with the results of plane geometry 
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e on the particle creation 



23] 



For the case I ^ 0, the resonances are shifted to noninteger values of the ratio w/uqi(0), 
since the eigenfrequencies uii s are no longer equidistant, as can be seen in Fig. 1. In this 
case Eq. (jSJ) has no analytical solution, so we can not write an obvious closed expression for 
the coefficients cf ss ,. In Fig. 2 we plot (r Q — r») \cf ss r\ as a function of the ratio r /?Y As we 
can see, the cases I = and I ^ exhibits the same behavior, i. e., as the distance between 
the shells decreases, \cf ss ,\ increases. In contrast, when I ^ 0, where c|/ M « ^ — (— l) s+s 'c^ ss ,- ) , 
the effective velocities v a — (— l) s+s do not exhibits an evident role in the amplitude of 
the coefficients cf, ss ,y This fact shows that the qualitative difference between the plane and 
spherical geometries appears essentially for I ^ 0, as evidenced from Fig. 1. 

Under the law of motion (J22J), we analyze the average number of particle creation in the 
four different cases mentioned in the introduction: when (a) only the inner shell oscillates 
(ej = e and e Q = 0); (b) only the outer shell oscillates (e« = and e Q = e); (c) both shells 
oscillate in phase (e« = e Q = e); and (d) both shells oscillate out of phase (e$ = — e = e). 

In Fig. 3 we present the plot of the ratio Ni ms (t) / (ewt) 2 against w/ujqi(0), under the 
resonance condition {w = lui ss '), for all four cases and for few values of I and s. We find that 
the principal resonance — which maximizes Ni ms (t) — occurs when w = 2su>i ss (0) (s' = s) 
in cases (a), (b), an (c), as expected. However, in case (d), this resonance can be shifted for 
the value s' = s + 1, depending on the ratio r /rj which, for I = reads 



This result follows directly from Eq. (|26|) . Note that the case where only the outer shell 
oscillates produces a larger number of particles than that where only the inner one oscillates. 
This can be observed directly from Eq. (1251) under the assumption of only one oscillating 
shell, where the rate Ni m8 (t) is proportional to the velocity of the moving shell. We also note 
from Fig. 3 that as I increases, the particle creation rate in this mode decreases accordingly. 
This result is expected, since the energy of a given mode increases with I. By its turn, for 
the cases (c) and (d) we can have a larger or a smaller number of particle created in the 
cavity, depending on the chosen resonance: if s + s' is an even number, the case (c) will 
present a smaller number of particles than that all other cases while the case (d) will present 
the larger number of particles among the four cases. If s + s' is an odd number, the opposite 
situation will occur. 

In Eq. ( 125]) we also observe that no particle can be created, even under the resonance 
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condition, under the following constraint 




> 1, 



which for I = reduces to 




s+s' 



>1, 



showing that, if s + s' is an even (odd) number and €i/e Q > (< 0) this condition is not 
satisfied and we will always have particles created in the cavity. 

We stress that since the expression for the average number of particle creation was ob- 
tained in the frameworks of the perturbation theory, for resonant breathing modes it is valid 
only in the short-time approximation, e a zut <C 1. Therefore, the analysis performed above 
cannot predict the real number of the created particles in the long time limit. 

V. CONCLUDING REMARKS 

Here we have considered two concentric spherical shells that are allowed to move, and 
we analyzed the particle creation within the DCE for the massless scalar field confined in 
the cavity. The particle creation were computed for an arbitrary law of radial motion of 
the spherical shells, using two distinct methods: the density operator of the system and 
the Bogoliubov coefficients. We applied our general results to the case of an oscillatory 
radial motion of the spherical shells, associated with breathing modes, identifying the reso- 
nance conditions were the number of particle creation is more significant. Analyzing these 
resonances, we have noted that qualitative differences between the plane and spherical ge- 
ometries arise when I ^ 0. We have considered four distinct cases of the breathing modes: 
when only (a) the inner or (b) the outer shell oscillates, or both shells oscillate (c) in phase 
or (d) out of phase. As already emphasized, our resonant results are restricted to the short- 
time approximation e a tot <C 1. In conclusion,we believe that the present work is enlarging 
perspectives in the subject of the DCE, which is getting the interest of both,theoreticians 
and experimentalists. 
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Figure captions 

Fig. 1 (color online). Map of the solutions of the transcendental equation (III]). The 
colors correspond to different values of the number I: the black lines are for I = 0, the red 
ones for / = 1, and for Z = 2 the blue lines. The solid, dashed, and dotted lines correspond 
to the cases s = 1, s = 2, and s = 3, respectively. 

Fig. 2 (color online). Plot of (r D — r^) \cf ls , | against the ratio r /Yj. The solid and dashed 
lines correspond to s' = 1 and s' = 2, respectively. The black line correspond to / = and 
a — i, o. The blue and red lines are for a = i and a — o, respectively, both with 1=1. 

Fig. 3 Plot of Ni ms (t) I (ewt) against zu/u 01 (0), in the resonance condition, for the cases 
(a), (b), (c), and (d) for few values of I and s. We have set r Q = 2rj. 
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